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We have calculated the Bardeen-Stephen contribution to the vortex viscosity for uniaxial 
anisotropic superconductors within the time-dependent Ginzburg-Landau (TDGL) theory. We fo- 
cus our attention on superconductors with a mismatch of anisotropy of normal and superconducting 
characteristics. Exact asymptotics for the Bardeen-Stephen contribution have been derived in two 
hmits: (i) Isab < ^ab, Iec < ?c and (u) Iec > Cc ^Eab < Cab, where Isab, Iec and Cat, Cc are the 
electric field penetration lengths and the coherence lengths in the ab plane and in the direction of the 
c axis. Also we suggest a variational procedure which allows us to calculate the vortex viscosity for 
superconductors with arbitrary parameters ^ and Ie- The approximate analytical result is compared 
with numerical calculations. Finally, using a generalized TDGL theory, we prove that the viscosity 
anisotropy and, thus, the fiux-flow conductivity anisotropy may depend on temperature. 
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I. INTRODUCTION 

The existence of a non-zero electrical resistivity in 
type-II superconductors in the mixed state is connected 
with the motion of magnetic flux vortices. It can be ob- 
served in the presence of a sufflciently large transport 
current so that pinning is suppressed. In the stationary 
flux-flow regime the Lorentz force acting on an isolated 
vortex is balanced by the intrinsic viscous drag force: 



(jt 



V^i 



(1) 



Here ^o is the flux quantum, jtr is the transport current 
density, n is the unit vector along the magnetic field, 
V^ is the vortex velocity and 77 is a viscosity coefficient. 
As vortices move, the magnetic field in the sample be- 
comes nonstationary and a macroscopic electrical field E 
is induced, which is connected with the transport current 
via Ohm's law: E = jtr/c. For weak average magnetic 
fields, B ^ Hc2, where i?c2 is the upper critical field, the 
flux-flow conductivity a is 
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(2) 



The presence of a flnite conductivity implies that cur- 
rent flow is accompanied by dissipation. It has been 
shown that there are two main mechanisms of dissipa- 
tion: losses due to relaxation of the order parametepU 
and ohmic losses associated with normal currents flow- 
ing through the vortex core.*^ 

For an anisotropic superconductor Eq. (pi) is general- 
ized as follows: 
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with the z-axis along the magnetic field. The peculiar 
structure of the conductivity tensor is explained by the 



fact that the x-component of the electric field depends on 
the y-component of the vortex velocity, and vice versa. 

It can be seen from Eq. ((sl) that the flux-flow conduc- 
tivity is determined by the magnetic field and the viscos- 
ity tensor r}. A rigorous approach to the problem of vis- 
cosity evaluation has been first suggested by SchmicP and 
was later developed by Gor'kov and KopniiP (see also 
Ref. inifor review). Their method is based on the time- 
dependent Ginzburg-Landau (TDGL) theory. Within 
this model the flux-flow conductivity has been evaluated 
for isotropic superconductors in several papers.!^"^. Both 
viscosity components due to order parameter relaxation 
and ohmic losses (frequently called the Bardeen-Stephen 
contribution) have been derived. 

Theoretical studies of free flux flow in anisotropic ma- 
terials have been stimulated by the discovery of high- 
temperature superconductors which appeared to pos- 
sess rather strong anisotropy. A number of papers 
have addressed this problem using different models and 
approximations.™3JJ xhe procedure of viscosity calcula- 
tion can be essentially simplified in the limit of dirty uni- 
axial superconductors with the ratio sq = mc'^c/'m'ab'^ab 
equal to unity. Here ac, Tab a-nd mc, niab are the nor- 
mal conductivities and Cooper pair masses in the direc- 
tion of the anisotropy axis c and in the perpendicular 
afe-plane, respectively. The condition sq = 1 allows to 
reduce the problem of anisotropic vortex dynamics to 
an isotropic one by means of a scaling transformation.Ml 
Yet this is not true in the case sq ^ 1, i.e, for a mis- 
match of anisotropics of Cooper pair masses and normal 
conductivities. Such a mismatch is theoretically possi- 
ble in the relatively clean limilP2 and it may have been 
experimentally observed in a new class of Fe-based pnic- 
tide superconductors. According to Ref. il3i the ration 
Oabl '^c in PrFeAsOo.7 is close to 120, whereas mdmab in 
the same compound is about 25, as determined in Ref. 
[14] from upper critical field measurements. In Refs. [TSl 
and [TBI anisotropics of the same order in Bai_xKxFe2As2 
are reported. However, existing experimental data for the 



pnictides are contradictory. In Ref. [T7]a relatively low re- 
sistivity aniso tropy in BaFe2As2 is given: Oab/cfc ~ 2 — 3. 
In some workaSUSl g^^ anisotropy mismatch has not been 
clearly detected. Previous calculations of the vortex vis- 
cosity tensor accounted for the anisotropy mismatch only 
on the basis of a simplified model of a step-like order pa- 
rameter profile within the vortex core W Of course, a de- 
tailed comparison with experimental data demands these 
calculations to be generalized for a more realistic order 
parameter profile. 

In this paper we evaluate analytically the viscosity ten- 
sor for a realistic gap profile within the core, focusing our 
attention on the case sq 7^ 1 and considering both stan- 
dard TDGL model and its ge neralization for supercon- 
ductors with a finite gap.li^ESlln section |ll] we derive the 
basic equations foll owin g the approach of Gor'kov and 
Kopnin.Ei In Section 



III 



we develop approximate methods 
based on different assumptions about the ratio of the elec- 
tric field penetration dept h to the coherence length. The 
results of preceding workaSZEl are revised and improved. 
In the end of this section we consider a variational prin- 
ciple which provides us with a simple general relation 
for the Bardeen-Stephen contribution. In section [TV] the 
problem is considered in the framework of a generalized 
TDGL theory. Here we derive our main result: we predict 
that the flux-flow conductivity anisotropy may depend 
on temperature in superconductors with the parameter 



II. BASIC EQUATIONS 

Following Gor'kov and Kopnin,^ we start the analysis 
of vortex motion with the TDGL equation for the super- 
conducting order parameter tp: 
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Here F is the usual GL free energy, 7 is a relaxation 
constant and A and $ are the vector and scalar poten- 
tials, respectively. We consider uniaxial anisotropic su- 
perconductors, so the Cooper pair mass to is a tensor 
with components rriij = mabi^ij + M^'i^'j) where u is the 
unit vector along the c-axis, ^ = rric/mab — 1- Eq. Q is 
supplemented by the equation for the current density 
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FIG. 1. The coordinate frame 



Here 9 — aig{ip) and (7„ is the normal-state conductivity 
tensor with components Unij — Oab^ij + [oc — (Tab)vii'j- 
For simplicity we will consider only superconductors with 
a large Ginzburg-Landau parameter n = A/^ 3> 1, where 
A is the London penetration length and ^ is the su- 
perconducting coherence length. One can prove that 
|2eA/c| / \hV9\ ^ 1 at distances much smaller than A 
from the vortex axis in the gauge where A = on the 
vortex axis and divA = 0. Imposing the additional con- 
dition l%/\£, ^ 1, where Ie is the electric field penetra- 
tion depth (see Eq. (12)), one can neglect the term 
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in Eq. (|6|. 

Let us consider the orientation of the internal mag- 
netic field at an angle Lp to the crystallographic c-axis. 
We choose the coordinate frame (x, y, z) with the z-axis 
coinciding with the vortex axis and with the c-axis lying 
in the a;z-plane (see Fig. nl). In this frame the functions 
■0 and $ do not depent on z, and the tensor i) is diagonal. 

The derivation of the force balance equation (fTl) and 
the explicit expression for the viscosity may be found in 
Refs. |S] and [HI However, in Appendix |^ we outline the 
calculations for the reader's convenience. 

The components of the Bardeen-Stephen contribution, 
7}oh, are given by 
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divj = 0, 
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where 
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rix = [TO((/3)/TOob]^^^(r?oh)xa: 
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dxidyi. (8) 
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(11) 

^(j6 and Isab are the coherence length and the electric 
field penetration depth in the a6-plane, respectively: 
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(12) 
The function f{p) describes the profile of the dimension- 
less order parameter modulus in a static isotropic vortex. 
This function satisfies the relation 



pdp \ dp J p^ 



(13) 



which follows from Eqs. (A8| and (A17). The boundary 



conditions are /(O) — 0, /(oo) = 1. The functions $2^ 
and ^y in Eqs. (l7| and (IS]) should be determined from 
the linear equations 
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where 



Siifi) = 1 + 
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1 + p cos^ (^ 



>0. 



(14) 
(15) 

(16) 



The electric potential can be expressed in terms of $2 



and $j, via 



This limiting case is more close to gapless superconduc- 
tors with a high concentration of magnetic impurities 
which are characterized by the ration ^/Ie — \/l2. The 
conditions (19) impose the following restrictions on the 
parameters s and u: u ^ 1, s ^ w^. In this subsection 
we will analyse the case s < 1. The case 1 ^ s ^ u^ can 
be considered in a similar way by dividing Eqs. ( 14 1 and 



( 15 1 by s. We shall search the asymptotics of the viscos- 
ity when u — >■ 00 neglecting small terms of order higher 
than w~2 (however, it will be shown that one should keep 
terms of the order of m"^). 

Our approximation is based on the fact that the char- 
acteristic length scale for the functions ^^ and ^y is m~^. 
Hence, the unknown functions reach their asymptotic be- 
havior at distances p <S^ 1 from the vortex axis, where the 
order parameter profile f{p) is well approximated by the 
first several terms of its Taylor series: 

f{p)-k2p''+hp'' + hp''. 



We substitute this expansion into Eq. ( 14 ) and introduce 
the new variables p = p^s/u, ^^ — ^^^^Fi"^: 



s- 
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k2P^ + ki- 



$. - -| ) . (20) 



Further the tilde over x and y will be omitted. The solu- 
tion of Eq. (20 1 can be expanded in the powers of u"^: 



$, = $i°)-fu-i4i)-)-i?„ (21) 

where $i and <^x satisfy the following relations: 

s ^ „ = kap'^^l-' - k2y, (22) 
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dx^ 



$ = $. 



''^V.x-'^yVjp'j^. (17) 
rUab / baab V 2TOab 



Note that there is a relation connecting the compo- 
nents ?7a; and ijy-. 



r]yis,u) =r]x 



1 u 



In the next section we consider some limiting cases. 
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= fc2P^$W+fc4pM'J'i°'-4). (23) 



and Rx is a remainder term. It is proved in Appendix [B| 
that an analogous expansion can be made in the integral 
in the rhs of Eq. ([t]): 



Vx = -2—-jh 
b 
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(24) 



III. APPROXIMATE SOLUTIONS 

A. The Zb < ^ limit 



where 



hx{s) ^ j k^y U^^^ - ^\ dxdy, 



(25) 



Consider such materials that the electric field penetra- 
tion length is much smaller than the coherence length: l2x{s) 



^Eab ^ Cab J ^Ec ^ Cc 



(19) 
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(26) 



The viscosity component ijy can be calculated similarly 
'/i,(.s) L 
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(27) 



Using Eq. ( 18 1 we obtain 



hy{s) = /l.(s-l)Vs, hy{s) - hx{s-^)s. (28) 

In principle, the functions Iix{s) and l2x{s) can be de- 
termined by numerical calculations, however, in section 
|III C| we present some analytical expressions for these 
functions. 

In Ref. 9 the u ^ 1 limit was considered using the 
Bardeen-Stephen model. ■^ This approach is essentially 
based on the assumption about a step-like order param- 
eter profile within the core and does not allow to obtain 
a leading term of the order of u 

The particular case s 

number of works mentioned abovc^l^t^ It corresponds 
to isotropic superconductors, or anisotropic superconduc- 
tors with no anisotropy mismatch: {fricCTc)/ {'niab'^ab) — 1- 
If s = 1, Eqs. (22 1 and (23) can be solved exactly: 



^ in the expansion (24 1. 
1 has been considered in a 
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After some integration we obtain a simple relation for the 
viscosity: 



Vx =r]y = 2TT — jha2{u), 
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(29) 



(30) 



Here the value ^/k2 = 0.583 was taken from Ref. [71 and 
the relation ^4 = — fc2/4 follows from Eq. (13 1. 



It is appropriate to recall here the result obtained by 
HuPl 



a2 



Ko(fa) 
Su ■ Ki((5u) 



(31) 



where Kg and Ki are the modified Bessel functions of 
an imaginary argument and 5 is a fitting parameter. Eq. 
(31 1 was derived from the exact solution of Eq. (A12) 



with an approximate order parameter profile: 
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(32) 



According to Schmid'^ and Hu, the optimal value of 5 is 
v2 which follows from a variational principle. We can 
compare different values of a2{u). When u — \/l2 Eq. 
^ yields a2 = 0.186, Eq. ^ yields az = 0.158, while 
the numerical result is a2 = 0.159.'^ Our formula gives 
an error less than 1%. If we keep only the term of order 
u~^ in Eq. (30 1, we will get a 6% error which increases 



with decreasing u. 



B. The Iec > ic limit. 

Consider the range of parameters s ^ m^ and u > 1. 
In terms of ^^ , 1^ and Lp these conditions read 

^£c>Cc, lEab<^ab, COS ^ < ^^^ • 

^abl;ab 

Thus, the magnetic field must make a small angle with 
the a6-plane. 



When s 3> u^, the term u^^x in Eq. (14 1 is negligi- 
ble compared to yi/ p\ in the region p ^ \/s/u, so we 
immediately obtain from Eq. ([7]) 

r]x ^ Ins/u^. 

More complicated calculations, which can be found in 
Appendix [C] yield 

r^x = 2TThj^-^ (\n ^ - 1.475") , (33) 



^2^/17^ (In ^-0.475 
\ u 



(34) 



Note that in Ref. 9 in the u <ti 1 limit similar expressions 
containing Inw^^ have been derived. This similarity is 
not accidental: the presence of the logarithm ln{lE/£,) is 
a characteristic feature of the Ie ^ ^ limit. 



C. A variational principle 

In this subsection we suggest a simple variational pro- 
cedure for the calculation of the viscosity tensor. Ac- 
cording to Ref. [5] a general expression for the dissipation 
function W[^] reads: 



W[<^] = V$(T„V$ + 
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h^ + 2Je$V 
at 



(35) 



The electric potential should be found from Eq. (A12) 



which can be viewed as a condition of zero variational 
derivative of the functional 



Jz=0 



Thus, the minimum of the functional above equals the 
loss power per unit length of a moving vortex: 



VlvVl = min / W[^]d^p. 

* Jz=0 



(36) 



This relation allows us to apply the direct variational 
method to our problem. 



For the sake of convenience we rewrite Eq. ( 36 1 in 



our rescaled coordinate frame separately for both com- 
ponents of the Bardeen-Stephen contribution: 



r)x=fi{s,l,u), r]y^fi{l,s,u), 



(37) 



fi{sx,Sy,u) = 2'-f-jhu'^ ■ minj 
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(38) 



Eqs. (37 1 - (38 1 have two important consequences. First, 



the viscosity tensor is positively defined when Gab > 
and (Tc > 0. Second, the components of f)oh increase as 
the conductivity increases: 



dcFab 



>0, 






>0, 



x,y. 



We can obtain an upper estimate for the viscosity com- 
ponents if we substitute a trial function into Eq. ( 38 ). In 



order to find an appropriate trial function consider the 
exact equation for (p: 






u^cj)- 



y 



fip)- (39) 



The solution of this equation is an even function of x and 
an odd function of y, so its Fourier series has the form 



y^ (P2n+1 
n=0 



02ri+i(p)sin(2n + l)x, 



(40) 



where x is the polar angle in the xy plane. When p is 
sufficiently large, (f> « y/{v?p'^), that means, that the 
series in Eq. (40) contains only the first term. Thus, the 
trial function 



#(p) 



smx 



(41) 



has the correct parity and the correct asymptotics. Let 
us substitute this function into Eq. ([38|): 



2 — jniTU mm 
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where 
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The differential equation for (j) is 



^IMh^/^")i^'* 



0. (43) 



Note that we obtain exactly the same equation if we sub- 
stitute $a; = (j){p) sinx into Eq. (14) when s — 1 and 
u — u. This means that the trial function (41) reduces 



our problem to an isotropic one. Unfortunately, an ex- 
act solution of Eq. (43) is unknown. However, Schmid'^ 



found a solution with an approximate order parameter 
profile [see Eq. (|32|]: 



Ki{u6)5 - v/(52 + p2Ki(M^,52 + p2) 

(5Ki(u(5)u2p 



Using this function and the expression ( 32 ) for / we can 



calculate the rhs of Eq.(|42|: 






(44) 
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We take 5 = /'(0)~^ — k^'''^ in order to obtain the cor- 
rect asymptotics when u -^ oo, s = 1 (this asymptotics is 
determined by /'(O), see subsection III A). Finally, com- 



componcnts of ?)' 



bining ( 37 ) and ( |44[ ) we get approximate relations for the 
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No restrictions on the parameters s and u are implied 
here. 

Let us check if these relations are in accordance with 



the results from subsections III A and IIIB Expanding 
r]x in the form (45 ) in the powers of vT^ when u ^ 1 and 



s < 1 we obtain the following expressions for the coeffi- 
cients /la; and Iix, which were introduced in subsection 



III A (see Eq. (24)): 



When s = 1 



TTy/k^y/s + 3 



7r/c2(s + 3) 
hx[s) = ■ (47) 



Vx=Vy = 2An^^+Oiu-') 
b u 



which should be compared with Eq. (30). The perfect 



agreement between the exact and approximate result is 
not surprising, because the trial function ( 41 ) is the exact 



solution of our variational problem in the isotropic case. 
In order to check whether Eqs. (45) and (46) are ap- 



plicable for s ^ 1 we used numerical calculations. We 
solved Eq. (39) in the region a; > 0, y > with the 



boundary conditions 



dx 



0, 



x=0 



y=o 



A sufficiently large 450 x 450 mesh with a 0.03 x 0.03 unit 
cell has been used. The numerical algorithm applied was 
the method of steepest descent. After the determination 



of the function 4>{p) numerical integration has been per- 
formed. 



When s = Eqs. B7| and (28 1 give 



\a\ ^1.59 ^, 2, 
b u 



\a\ ^0.92 ^, 

r]y = 2^-^jh +0{u- 

b u 



These analytical expressions are in a good agreement 
with the asymptotics derived by numerical calculations: 



\a\ ^1.58 
u 



0{u-% 
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\a\ Um ^, n, 
u 



When s > u Eqs. (45 1 and (46) give 



r\x 



Vy 



b u 



which coincides with the main logarithmic term in Eqs. 
(Isl and (IMl). 



One can see that the agreement between the exact and 
approximate asymptotics is quite well. This is a strong 
argument in favor of the applicability of Eqs. (451 and 
(461 for intermediate values of s and u. 



In Fig. [2]we plot the analytical and numerical Lp depen- 
dencies of the diagonal components of the full viscosity 
{fj^VpO + 'noh)- 



IV. TEMPERATURE DEPENDENCE OF THE 
VISCOSITY ANISOTROPY. 
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FIG. 2. The ip dependencies of the viscosity components. 
Solid lines correspond to analytical results (Eqs. ( |45[ ) and 
(46l), dashed lines show the results of numerical simulations. 
Here rixx and r\yy are measured in the units 770 = ^l\o\ /h, 

iab/lEab ~ V 12. 



Within the framework of the TDGL equation (B the 
viscosity anisotropy rj^x/Vyy does not depend on temper- 
ature. However, the region of applicability of Eq. ^ is 
limited by gapless superconductivity. In this section we 
consider a more general approach based on the general- 
ized TDGL equationJiSl (see also Ref. [T^lfor review): 
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Here T^ is the critical temperature and Tph is the electron- 
phonon mean free time. In the isotropic case Eqs. (48) 



and (49) are vahd for dirty superconductors, when the 



temperature is close to Tc and variations of the order 
parameter in space and in time are sufficiently slow. 

The main relations for the viscosity can be derived in 
same way as described in section [IT] As a result, we find 
that the viscosity still comprises two terms representing 
two mechanisms of dissipation, but the viscosity compo- 
nents undergo some changes. For example, Eq. (A19) is 
modified as follows: 



iai„ ruf^^ 



(ry^o),, = 2nhj'-^6,, / ( ^ ) p^l + qp{p)dp. (50) 



In order to obtain the counterparts of Eqs. (l7|, (|8|, (14) 
and ( 15 ) one should make the following substitutions in 
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FIG. 3. Schematic temperature dependence of the viscosity 
anisotropy for tp — -k 12. The parameter e is of the order of 
unity. 



these equations: 
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i^x,y. 



(51) 



It can be seen from Eq. ( |49| ) that the electric field pene- 
tration depth is increased by a factor (1 + qY^'^ as com- 
pared to Eq. (12]). It may seem that at low temperatures 
we would reach the Ie ^ ^ limit, which has been analysed 
in Ref. [HI However, this is not quite true because of the 
different relative impacts of the two mentioned mecha- 
nisms of dissipation in the simple and generalized TDGL 
models. Within the simple TDGL theory the Bardeen- 
Stephen contribution and the relaxational term are of 
the same order of magnitude in the Ie ^ S, limit. On the 
contrary, in the generalized model the viscosity is domi- 
nated by the relaxational term at low temperatures (see 
below). 

It is obvious that all main relations from sec. IIIII can 
be derived again within the generalized TDGL theory, 
but they are slightly modified. For example, Eq. (29) 
now reads 



Vx ^ Vy ~ ^nh'y- 



u 



l + 2q 
8u2 



(52) 



Now consider the temperature dependence of the vis- 
cosity. The quantity q depends on the temperature T, 
and gr^ < 0. Hence, 



d ir]po)a 



d ir]po)yv 



<0, 



where rjo 



dT 770 dT 770 

h'^\a\/h. On the other hand, the modified 



Eq. ( 38 1 can be written in the form 
fj{sx,Sy,u) = 2'-^^hu'^mmJ 



^ \dxj 

dxdy, (53) 



if we leave Eqs. (37) unchanged. Hence, 

d ir]oh)xx ^ „. d {rjoh)yy 



dT 770 



>0; 



dT 770 



>0. 



At sufficiently low temperatures, when g 3> 1, s ^ 1 and 
M ~ 1, it may happen that 



{r]po)xx > (»7oh) 



{ripQ)yy > ('7oh)i 



Then the viscosity anisotropy is determined by the relax- 
ational term: 



Vx 



{Vpo)xx 



l + ^J. 



Vyy i'npQ)yy I + fJ.COs'^ if' 
Note that when s 7^ 1 772, 7^ ?7j^, so 



Vyy 



^ 



l + M 

1 + /i COS^ ip 



when g < 1. We have proved that within the generalized 
TDGL theory the viscosity anisotropy and the flux-flow 
conductivity anisotropy do depend on temperature. The 
schematic T dependence of the ratio rj^x/Vyy is plotted 
in Fig. |3] 



V. CONCLUSION 

By solving the time-dependent Ginzburg-Landau equa- 
tion we analyzed the viscous flux-flow in anisotropic su- 
perconductors. The Bardeen-Stephen contribution to the 
viscosity tensor 77 has been calculated in the /^ <C C ^md 
Iec ^ £,c limits. We emphasize that in these calculations 
we did not use any simplifying assumptions concerning 
the shape of the order parameter in a static vortex. We 
suggested a variational procedure, which allowed us to 
derive the relations (45) and (46) suitable for arbitrary 
electric field penetration lengths (lEab and Iec): coher- 
ence lengths {S,ab and S,c) and orientation of the mag- 
netic field. Our results may be useful for interpretation 



of experimental data on flux flow conductivity in isotropic 
and anisotropic superconductors in weak magnetic fields 
{B « H,2). 

Viscous flux-flow has also been examined within a gen- 
eralized TDGL theory. We found that the viscosity 
anisotropy may depend on temperature and, thus, the 
flux-flow conductivity anisotropy may be altered by heat- 
ing or cooling the sample. We hope that this effect will 
be observed experimentally in the near future. 
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APPENDIX A 



In this appendix we derive Eqs. ([t]), (|8|, (14 1 and 
(15 1. In Eq. (|4|) it is convenient to make a scaling of 

y, z' = z. We 



the variables: a?~= x{m{(p) /mab)^''^ , v' 
rewrite Eq. (HI) in the form 

_d\ip\ h^ 



7/1- 



dt 



2m„ 



TlVlM^I^+Se* 



V^\^\-m{V0r]-a\^P\-b\^\\ 

(Al) 

-V'd^l'v'e). (A2) 



2mab 



According to Eq. ([5|, the two-dimensional current j' 



{jx[m{if)lma 



satisfies the relation 
div'j' = 0. 



(A3) 



It follows from Eqs. (A2) and (A3) that 



71 



2e$ 



-V (a„V'$) 



dt ' J 4e 

where we introduced the tensor (7„ with components 



(A4) 



fJ. 



nx'y' 



^7iy' x' 



0, 



(A5) 



with s{lp) given by Eq. (16). For a moving vortex one 



should search the solution of Eqs. (Al), (A2) and (A3) in 



the form 1/' = ipip'-VLt), $ = <i>(p'- V^t), where Vl == 
{VLx['m'{(p)/mab]^^'^, VLy) and p' = {x',y'). We expand 
1^1 and in powers of Vl up to the first order term, 
assuming the vortex velocity to be sufficiently small: 

|^|«^o(p'-Vit) + 7^i(p'-Vit), (A6) 

0«0o(p'-Vii) + ^i(p'-Vit). (A7) 

Here ipo{p') and Oo{p') correspond to a static vortex. The 
functions V'l, ^1 and $ are of the order Vl- We substitute 
(lAel) and (TAtI) into Eqs. (lAll), (lASl) and (TaI: 



aiJo - b^l + -^ [V^ijo - M^'Oof] - 0, (A8) 



div'j^ = 0, 



Jo 



2ehtp^ 

rUab 



V9n 



(AlO) 






div'j'i = 0, 



a„V'$, (All) 



flOab 



s{ip) 



92$ 92^ 



dx'"^ dy''^ 



-1^1 (2e$ - Wl • V'0o 



(A12) 

Now we introduce some new notations: V'd = (dV')'0Oj 
9d = (dV')6'o, j^ — (dV')Jo, where d is an arbitrary 
vector. A simple equation connecting ^^ and 9d can be 
obtained by applying the operator dV to Eq. (A8): 



- ailJd - ih^l^d = (A13) 

The vector jjj satisfies the obvious relation div'jj; = 0. 
Let us multiply Eq. (|A9|) by -0^, subtract Eq. (|A13|) 



multiplied by "01 and integrate the resulting equation over 
a large volume containing the whole vortex. After some 
simple algebra and integration by parts we obtain 

-7;i/(ViV')^o0dd'r' 

= s / [(J'l + a„ v'$)v'0rf - I'^ve^] Sv' 

= I /(a„V'$)V'0,d3r' + A S^{]'^ed~i',e,)dS. (A14) 

Here S* is a surface far from the vortex axis. At large 
distances p' ^ S,ab we have 



Ji 



2eh\a\ 
hruab 



V'^i^jir, 



bmgb 
2eh\a\ 



(Jtr • P') + const. 



where j^^ is the transport current which is constant. If 
we calculate the surface integral in Eq. (A14| and make 



some sirnple transformations, we obtain the force balance 
equatiorP 



irh 



[d. (j;, X n)] = -27rjhid- Vl)J,^ {"tY Pdp 



2e 



I si^) 



dx'2 



dy'^ 



{d-Veo)<Pp', (A15) 



where n is a unit vector along the magnetic field. If 
we compare Eqs. (fTl) and (A15), we can see that the 
viscosity tensor in the frame (x', y' , z') should be defined 
as follows: 



N 
h 

where we introduced the function 



d.77'Vi = 2^7S(d.Vi)' 
92$ 



df 



pdp 



hcFgb 

^2e 



<^)-. 



dp ^ 
(d.V%)d2p',(Al6) 



2ma 



[V'20i - M^'Oo? - 2^-0 V'00 • V'^i] 
- a0i - 35^2^1 = -jh{VL^')^o 



(A9) 



f{p) = i^o{pLb)] 



(A17) 



9 



The components of the viscosity t) in the frame (x, y, z) 
are given by 



order of (wpo) ^ <C w ^, the integrals of higher-order 
terms are also negligibly small, hence 



(A18) 
The rhs of Eq. ( |A16 ) contains two terms, representing 
two mechanisms of dissipation. The viscosity due to re- 
laxation of the order parameter igSEl 



11x2 



J' - 



P>Po 



b U 






(B5) 



dx dy. 



{v'po)i3 =2Trhj—aiSij, 



fA19) ^^^ ^"^^ consider the component 773,1. In the new variables 
introduced in subsection III A Eq. ( B2 1 reads 



ai 



pdp = 0.279. 



The second term in the rhs of Eq. (AI6) defines the 



ohmic viscosity tensor 77^,^, which is to be evaluated: 



d • €h^L 



ha„ 



2e 



s{ip) 



dx'"^ dy''^ 



(dV'6'o)d>' 



'2„' 



Now, if we substitute $ in the form (17) into Eqs. (A12) 



(A20) 



and (A20) and switch to the coordinates (£1,2/1) (see Eq. 
([n])) we obtain Eqs. ([7|, ^, ^ and ([l5|. 



APPENDIX B 



In this Appendix we will derive Eq. (24). First, we 
divide the integral in Eq. (It]) into two parts: 



rix = Vxl +11x2, 



(Bl) 



a 



2.„ ^2/l 1^2^ 



Vxl = -2^7/1/ / (Pl; 2 

" Jpi<pa/V^i Pi 



Pi 



-^ ]dxidyi, 



(B2) 



Vxl 



..^.» 



f 



P<PO 



y 



$. 



y 



dx dy. 



Now we estimate the term R^ introduced in Eq. ([21j) . It 
satisfies the following relation: 



(B6) 



yfl^ I „ o fix 

dy^ +« 



dx2 



^f' (^) ^^ 



PU,)-k: 



\V^ 



«2P 



-f {^u 



1 4 

K4P 



k2p' 



^^^j. 



(B7) 



Note that when p <^ ^Ju the source in the rhs of ( B7 ) can 
be presented as S{x,y)u~'^, where S{x,y) is some func- 
tion independent of u. Since (B7) is a screening equation. 



the function Rx{x,y,u) for small p does not depend on 
the behavior of the source in the area of big p and can 
be presented as Rx ~ Rx{x,y)u^^. On the other hand. 



when p '3> I the derivatives in the Ihs of Eq. (B7) are 
small, hence in the area 1 ^ p <?C ^/u 



Rx 



1 



kep'^ 



$(0) 



y 



kip" 



$ 



(1) 



\Rx\ < 



const 



pu^ 



(B8) 



11x2 



-2^7^ 



fipi) 



2/1 



Pi>Po/\/u Pi 



U'^:r 



yi 
pI 



dxidyi 



(B3) 

where po = u^/^+\ 5 < 1/6. Note that the Ihs of Eq. 
(14) is small when pi > po/^/u ^ u^^/^, so it can be 
accounted for by perturbation theory: 



$r = 



yi 
u^pI 



+ *xl+*x2 + -, 



(B4) 



Now we substitute $3, in the form (21) into (B6) 



\ u u 



J' 

^2x 



'3x 



where 



/;,= / k2y{^^^^~^]dxdy, 



' P<PQ 



(B9) 



(BIO) 



^xl 



$x2 = 



92 d'^ \ 2/1 

+ S-, 



u'^PiPi) \9yj dxlJ pI' 



92 



92 



P{pi) \dyi dxl 



yi 
2 ' 
Pi 



The main contribution to the integral in Eq. (B3) is 



determined by small pi. The integral of '^x2 is of the 



I'2x 



y_ 

P<PaP 



k,p^(<,io,_y_^yk,pHi^^ 



dx dy, 
(Bll) 



k2p __ k^p 

U li^ 



k^Pl 



^3x -jp<pa[f [^) 



K-J^) 



*r-S)Srfa;d2/ 



-^^^^dxdy 



(B12) 
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One can easily prove that 



^ P^ k2P^ \dy^ dx^ J p^ 



^' siH+0(p-^), (B13) 



k,p'U(^^-y-]+k,pHi'^^oip- 



(B14) 



From Eqs. (|B13f , JB14|) and (lB8|) we can see that all 
integrals in Eq. (B12| are of the order p^/u^- Thus 



\Izx\ ^ u ^, S O it c an be neglected. Also we can integrate 
in Eqs. (BIO) and (Bill over the whole xy plane, since 






(B15) 
(B16) 



Finally, taking into account Eqs. (Bl), (B5) and (B9) we 



obtain Eq. (24) 



APPENDIX C 



In this Appendix we consider in detail the derivation of 



Eqs. (33 1 and (34). We will present here the calculations 



for the rjy component, since the calculations for the rjx 
component are less complicated. First, we rewrite Eq. 
( 15 ) in the form 






u'^ + fiPip) 






u'%{l-f{p)). (CI) 
The index "1" is omitted. It will be proved below that the 



terms in the rhs of Eq. ( CI ) give a small contribution to 



the viscosity, so they can be neglected. Then the solution 
of Eq. ( CI I has the form 



$. 



%o = 



-oo / f2l „l 



x'f\p')cM-u\x~x'\l^-s) 



dx'. 



(C2) 

where p' = \/x'^ + y^ . Consider a quantity j/o in the 
range 1 <C j/o *^ \fs/u (for example, j/o = s^/^/v}/"^). 
We divide the integral in Eq. (Is]) into three parts: 



riy 



2^7/^ 



4l<yo'^'(/');^"'*^^^^y" 



Using the inequality 



f{p) < 



P 



(C4) 



where Ai is some constant, we can estimate the first in- 
tegral: 



\y\<yQ 



f{p)^u^%dxdy 



^ .you /, y/s 
< const — ^ In 



<1 

(C5) 



Here and further "const" denotes a constant independent 



of any parameters. The second term in Eq. ( C3 ) has the 
following asymptotics when j/q ^ 1- 



f {p)—dxdyfiiTrlnya + Cy. (C6) 

\y\<yo P 



The constant CL will be evaluated below. The third inte- 



gral in Eq. ( C3 ) can be simplified if we take into account 
that p > yoTp > J/o and j/o ^ Ij so we can substitute 
unity instead of p : 

Iy>yof'ip)^{^''^y-^)d^'iy 

'^ Jvn yS i-oo .t2+4„2 "2; — 1^^ \dy 



^yo \ \/s J — oo x^+Ay^ 

4y 



2y 



f /o°°'^^/r'^y( 



sx^ /u^-\-4:y^ y 



«^ -f r 1^ if^e""'^^ = f 1^2/0 - f In 4^ + f C, (C7) 
where C is the Euler constant: 



In a; • e^'^dx w 0.577. 



Using (C3) - (C7| we obtain 



riy = 27rfi,7 



TT 



b \2"' Au^ 
The component rjx can be calculated in a similar way: 



(C8) 






a 



- C , (C9) 



C, = lim / f{pf-^dxdy-Ti\nyA. (CIO) 

Now we evaluate Cx and Cy. Here the constant g'^ 
from Ref. 7 will be useful: 



9i = 



f{p)- 



l+p2 



p-^dp = -0.3982. (Cll) 



It is easy to check that 



Cx = ■^94+ lim / y'^p *dxdy ^ n (g'^ + ln2- -]. 

yo^°°J\y\<yo,p>yo V ^J 



Similarly, 



Cy^Tr[g'^ + \n2 



1 



Finally, the components of the viscosity take the form 
rjx - 27r;.7y ^^ V + ^^ " ^ " ' ^^^^^ 



11 



r,y = 2nhj^-^L^+g',-C+l]. (C13) 



4'^=Jfip)^u^R'^^dxdy. 



(C17) 



If we substitute C and g'^ with their numerical values, we 
obtain Eqs. (|33| and pt). 



Now it is necessary to prove our assumption concerning 
the rhs of Eq. (CI I. Consider it as a perturbation. The 
first order correction to the approximate solution <i>j,o has 
the form 

(I) 1 = i?(i) + i?(2) 



We will show that 



r(i) 



< 1 and 



r(2) 



^ 1 when s':^ V? . 



A simple estimate for |$yo| can be obtained with the 
help of (ICll): 



$yo| < const - 



In 



^ 



Ly/s 



Using the inequality 



P 

A2 + P 



2 < f{p) 



(C18) 



^(1)^/ 'u<^y,(.'^y)[l^-P{p')]^^^(^_u\.-yy^^^,^ 



y 2u^/s 



2^5 



1 /•+°° 92<I>,o , , 



(CM) 



dy"^ 



The contribution of <i>j,i to r]y is equal to 



(x',y)exp(-^^L^^d^'. 
(C15) 



where 



Aryy 



-2^7/.(X 



7-(2) 
1/ 



/(I) 






(C16) 



and Eq. (CI8I we can estimate ly 



(1). 



y 



< const ^ In ^^ < 1. 
Vs u 



For all x' and y we can write 



1 



d^ 



dy"^ \x'^ + j/2 



fip' 



< 



const 



(A3 + x'2 + y2)2 



whence 



r(2) 



< 



/ 



1 1 , u|:r-a;'| \ / u\x" -x'\ 
I X I CXp I ^ — y= 1 CXp I ^ 



^ const jjj \/£ r 
— s u J 

when u> 1. 



cxp 



a;"2+A, 



(A3+2;"2)(Ai+2:2) dx dx dx 

da;'dx" < £^ In ^ < 1 
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